Goos-Hanchen and Imbert-Fedorov shifts of polarized vortex beams 



Konstantin Y. Bliokh, 1 ' 2 * Ilya V. Shadrivov, 1 and Yuri S. Kivshar 1 

1 Nonlinear Physics Center, Research School of Physical Sciences and Engineering, 

Australian National University, Canberra ACT 0200, Australia 
Institute of Radio Astronomy, ^ Krasnoznamyonnaya Str., Kharkov 61002, Ukraine 
* Corresponding author: k.bliokh@gmail.com 



oo 
o 

O 

(N 

O 

Q 

m 

(N 



O 

q : 
o 

>v 

Oh! 



(N 
> 

in 
in 
oo 
(N 



OO 

o 

X 
S3 




Compiled December 23, 2008 
We study, analytically and numerically, reflection and transmission of an arbitrarily polarized vortex beam on 
an interface separating two dielectric media and derive general expressions for linear and angular Goos-Hanchen 
and Imbert-Fedorov shifts. We predict a novel vortex- induced Goos-Hanchen shift, and also reveal direct 
connection between the spin-induced angular shifts and the vortex-induced linear shifts. © 2008 Optical 
Society of America 

OCIS codes: 240.3695, 260.5430, 080.4865. 



Reflection and transmission of plain waves at an inter- 
face separating two homogenous isotropic media are de- 
scribed by the well-known Sncll's and Frcsncl equations. 
However, the problem becomes tricky for confined beams 
with a finite spectral distribution. Interference of partial 
plain waves propagating at slightly different angles and 
obeying individual Snell and Fresnel formulas, results in 
such effects as the longitudinal Goos-Hanchen (GH) [1,2] 
and transverse Imbert-Fedorov (IF) [3, 4] shifts, which 
displace the output beams within and across the propa- 
gation plane, respectively (see Fig. 1). 

While the GH shift was explained and calculated soon 
after its discovery [2], the transverse IF shift was asso- 
ciated by significant controversies over about 50 years. 
The direct calculation of the IF effect was first performed 
for the reflected beam [5-7] and later generalized to the 
transmitted beam [8,9]. It was shown that the IF shift 
is closely related to the spin angular momentum carried 
by a polarized beam and conservation of the total angu- 
lar momentum in the system [10,11] (see also [12-14]). 
Despite long history of the theoretical studies and experi- 
ments [4, 15,16] , analytical expression for the IF shift of a 
polarized Gaussian beam was derived in the correct form 
only recently [13,14], and these results have been con- 
firmed both experimentally [17] and theoretically [18]. 

In addition to the usual linear shifts, angular GH and 
IF shifts caused by the beam diffraction have been de- 
scribed [14,18]. Recently, this description of the IF ef- 
fect has been extended to the case of higher-order vortex 
beams carrying intrinsic orbital angular momentum [19]. 
The vortex-induced IF shift is proportional to the vor- 
tex charge, but it also significantly depends on the beam 
polarization. Such IF shift was calculated [19] and meas- 
ured experimentally [20] for p and s polarizations only, 
when the spin IF effect vanishes. 

The main purpose of this Letter is twofold. First, we 
derive explicit analytical expressions for both linear and 
angular IF and GH shifts in the most general case of 
an arbitrarily polarized vortex beam. We unveil a direct 
relation between the vortex-dependent IF shift and the 



Fig. 1. (Color online) Geometry of the beam reflection 
and transmission at an interface. Projections of the GH 
shifts, ~(X r ) cos 9, (X 1 ) cost?', and the IF shifts, (Y r '*), 
are shown by white arrows. 



angular GH shift [18] and predict novel vortex-induced 
GH shift related to the angular IF shift. Second, we ver- 
ify all theoretical results by direct numerical simulations. 

We consider the reflection and refraction of an opti- 
cal beam at an interface separating two media, as shown 
in Fig. 1. In addition to the coordinate system (x,y,z) 
attached to the interface, we employ the coordinate sys- 
tems of individual beams (X a , Y a , Z a ), where a = i,r,t 
denotes incident, reflected and transmitted beams, re- 
spectively. The Z a axis attached to the directions of the 
a-th beam as determined by the Snell's law. The inci- 
dent beam propagates in the (x, z) plane, so that Y a = y 
(see Fig. 1). We also define the wave number in the first 
medium, k, the angle of incidence, 0, the angle of refrac- 
tion, 0' = sin~ 1 (n _1 sin^), as well as the relative permit- 
tivity e, permeability fi, and refractive index n = ^/ejl 
of the second medium. 

We assume that the incident beam is a uniformly 
polarized Laguerre-Gaussian beam with the waist lo- 
cated at the interface, so that the transverse (X 1 , Y 4 )- 
componcnt of its electric field has the form 
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Here £j n ,m is the Laguerre-Gaussian solution of the scalar 
parabolic wave equation which contains an optical vor- 
tex of charge m: £„ jTO cx (X 1 + i sign (m) Y l ) '™ l [21], 
whereas polarization components form the normalized 
Jones vector in the basis of p and s modes: (e||,ejJ T . 
N| 2 + |eJ 2 = l. 

The main characteristics of the reflected and trans- 
mitted beams are determined by the Fresnel coefficients, 
i?H, R±, Tii, and Tj_. The amplitude reflection and trans- 
mission coefficients i?, T and the corresponding energy 
reflection and transmission coefficients Q T ' 1 are given by 



Here we denote and as the X a and Y a wave 
vector components in the a-th beam, D = kw 2 /2 is the 
Rayleigh length, wq is the minimum beam waist, and 
the derivatives of the energy reflection and transmission 
coefficients (7) are explicitly calculated (for the p and s 
polarizations) in [19]. 

Now we consider a polarized incident beam with vor- 
tex, m ^ 0. Such beam carries intrinsic orbital angular 
momentum (AM) L 4 ~ mZ 1 (in units of fi) per one pho- 
ton [21]. From the Snell's laws it follows that the AMs 
of the reflected and refracted beams are given by [19,22] 



R = ^\R VA] \ 2 + \R±.e±_\\ T = ^\T ]n \ 2 + \T x e±\*, 

Q r = \R\ 2 , Q^^dTf, Q r + 4 = 1, (2) 
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The Jones vectors of the secondary beams are 
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To derive the lateral shifts of an arbitrarily polarized 
vortex beam, we use the basic results obtained for polar- 
ized Gaussian beams (m = 0) . In the regime of the total 
internal reflection, sin# > n, the Fresnel reflection coeffi- 
cient become complex: R\\,± = exp(i<^ii and the GH, 
(X), and IF, (Y), shifts are described by the Artmann 
and Schillings formulas [2,5]: 
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where a — 2Im(ejjej_) is the helicity of the incident beam 
(degree of circular polarization), \ = 2Re(ejjej_) is the 
degree of linear polarization inclined at 7r/4 with respect 
to the incident plane, and S = ip± — <p\\ . We note, that 
the shifts here and below in the paper are given in the 
coordinate system of the respective beam. 

In the regime of partial reflection and transmission, 
sin (9 < n. the Fresnel coefficients are real, and the GH 
shift vanishes, (A" ,%t ) = 0, while the IF shifts arc given 
by the equations derived in [13, 14]: 
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In addition, the reflected and transmitted beams un- 
dergo angular shifts, which can be considered as the 
shifts in the wave vector space. Generalizing results of 
[14, 18], these shifts can be written as 
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As was shown by Fedoseyev [22] for a particular cases of 
p and s polarization, the IF shift of vortex beams con- 
sists of two contributions. The first one is a polarization- 
independent shift that comes from the difference in the 
z-components of intrinsic AM at the reflection and re- 
fraction which is compensated by the transverse shift 



producing an extrinsic AM [10-14]: (Y r 



(LI 



Ll)/ksm9. Taking into account that L l z = L l cos9, 
L r z = -L r cos9, and L* = L*cos0', one obtains [19,22]: 



(y)i = 0, (y*) 1 = ^tan0(l-n- 2 ). 



(10) 



The second contribution is essentially polarization- 
dependent. Here we demonstrate that it is directly re- 
lated to the angular GH shift (7). Indeed, as is shown 
in [18], the angular shift (7) induces an imaginary 
shift in the Gaussian envelope of the beam: X r ^ — > 
X r ' 1 - iD r /(K r /) /k (where D r x = D and L>* = 
D cos 2 9' I cos 2 9 [14]). It can be readily seen that this 
imaginary shift produces a real shift of the vortex in the 
orthogonal direction: 
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Here the coefficients j r — — 1 and 7* = cos 9/ cos 9' ac- 
count for the deformations of the vortex in the secondary 
beams: charge flip in the reflected beam and an elliptic 
deformation of the transmitted beam. As a result, the 
centers of gravity of the reflected and transmitted vor- 
tex beams experience the IF shift 

(Y r } 2 = mj(K r x )o, (F*) 2 = -m^^TOo. (11) 

The net IF shift of an arbitrarily polarized vortex 
beam is the sum of the contributions (6), (10), and (11): 



(Y r,t ) = (Y'.*) + (Y^x + (Y r >%. 



(12) 



This is the first main result that describes the total 
polarization- and vortex- dependent IF shift. 

Similarly to the vortex-induced IF effect (11) associ- 
ated with the angular GH shift (7), there exists a re- 
ciprocal effect of the vortex-induced GH shift caused by 
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the angular IF shift (8). In a manner, similar to the dis- 
cussion above, we obtain 

^9T {K y )o - (13) 



This is the second main result of this Letter. It predicts 
a completely new type of the GH shift, which occurs in 
the regime of partial reflection, in a sharp contrast to the 
usual GH effect (4). The vortex- induced GH effect van- 
ishes for p, s, and circular polarizations of the incident 
beam and reaches maximal values for linear polarizations 
inclined at it /A angles: x = ±1. 

It can be shown that the angular GH and IF shifts, 
Eqs. (7) and (8), acquire additional factor (1 + \m\) in 
the case of incident vortex beam: 



r,t \ 



l 



(14) 



The above equations (4)-(8) and (10)-(14) describe all 
the main shifts of polarized vortex beams with an axially- 
symmetric intensity profile. One may observe that the 
angular shifts (7), (8), and (14) fulfil conservation laws 
for the x and y components of the total linear momentum 
in the system [23]: 



f (K r x ) cos 6 + Q\Kl) cos 9' ■■ 

Q r ( K y) + Q l { K l) 
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whereas linear IF shifts (6), (10)— (12) fulfil conservation 
law for the z component of the total AM [10,11,13,14,22]: 

Q r (J r z - (Y r )ksm9) + Q t (J t z - (Y t )ksm9) = J*. (16) 

Here J z is the z-components of the total intrinsic AM 
of the a-th beam, which consist of orbital (vortex) and 
spin (polarization) contributions: 



ga ga = a a z a 



(17) 



where a 1 = a and a Til = 2Im(ejj' e r f) are the helicities 
of the corresponding beams. 

To verify our results, we perform numerical simula- 
tions of the problem based on the plain-wave decompo- 
sition of the incident polarized vortex beam and Frcsncl- 
Sncll's formulas applied to each plane-wave component 
of the beam spectrum. The numerically calculated co- 
ordinates of the centers of gravity of the reflected and 
transmitted beams in the wave-vector and coordinate 
space are compared with analytical expressions (6)-(8), 
(10)-(14) in Fig. 2, which shows an excellent agreement. 

In conclusion, we have calculated linear and angular 
GH and IF shifts arising upon the reflection and trans- 
mission of a polarized vortex beam at an interface be- 
tween two media. We have established a close relation 
between the spin-induced angular shifts and orthogonal 
vortex-induced linear shifts, and have predicted a new 
type of the GH shift for a vortex beam, which occurs in 
the partial-reflection regime when the GH shift for Gaus- 
sian beams vanishes. The value and sign of the vortex- 
induced shifts can be controlled by the vortex charge. 
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Fig. 2. (Color online) Numerical (curves) and theoreti- 
cal (symbols) results for the linear and angular IF and 
GH shifts of the reflected (r) and transmitted (t) vortex 
beam with complex polarization, as a function of the 
angle of incidence 9. Parameters are: n = 1.5, = 1, 
e 1 1 = l/v3) ej_ = (1 + i)/V3, and m = 1. For the 
wavelength of 632nm, (X r ) = 44nm for 9 = 40deg and 
(Y r ) = — 240nm for 9 ~ 90deg. These values are propor- 
tionally enhanced for higher m. 
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